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Physical Mechanism of the d — > d + zs Transition 
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We discuss the basic physical mechanism of the d — > d + is transition, which is the currently 
accepted explanation for the results of tunneling experiments into ab planes Using the first- 
order perturbation theory, we show that the zero-bias states drive the transition. We present various 
order-of-magnitude estimates and consistency checks that support this picture. 



It has now been firmly established |flj that the order parameter in cuprate superconductors has a 



the Andreev equation along each trajectory has a normalizable eigenstate of zero energy 



d-wave symmetry. An inhomogeneity may therefore scatter a quasiparticle between directions that 

o 

experience opposite signs of the order parameter. This effect is strongest for specular reflection off the 

' 110 surface. Here the sign of the order parameter along every quasiclassical trajectory changes sign 

> : 

upon reflection (see Fig. 1). As a consequence of this and the Atiyah-Patodi-Singer index theorem 
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These all contribute to the local density of states, and are the source of the zero-bias anomaly 

O ■ 

+^ \ observed in tunneling experiments [|§-0]. 

The tunneling data |]] show that the zero-bias peak splits as the temperature is lowered below 1 

1 : . , k „ .„„,,. . „-,,, . — „_,„.„,.. «,„„„, _ «,„ 

d — > d + is transition where a subdominant s-wave order parameter appears close to the surface 

> ■ 

^ , with a 7T/2 phase shift compared to the dominant ci-wave. Such a transition had been anticipated 
theoretically || prior to the experiment. Subsequent calculations of the surface phase diagram PJID|] 



included the effects of the surface roughness, and gave results that are in good agreement with the 
experimental data. 

These calculations are done within the Eilenberger approach to superconductivity fTI|| , which uses 



the quasiclassical approximation of the local Green function. This (so-called Eilenberger) function 



then satisfies a transport-like equation. Unfortunately, this approach is very formal, so it hides, 
rather than explains, the physics behind the calculations, especially when the equations are solved 



numerically in the imaginary-time domain fqjl^. As new experimental results |L3[ call the d + is 
interpretation into question, it is even more important to have a direct physical understanding of 
the transition. 



We believe the Andreev quasiclassical picture [14], which deals with single-particle states rather 



than local Green functions, contributes to such understanding by showing how the energy costs and 
benefits of the d — > d + is transition are distributed among the degrees of freedom in the system. As 
we mentioned above, due to the presence of the 110 surface, there are many zero-bias states (ZBSs) 
in the spectrum. The main purpose of this paper is to show that these states drive the transition. 
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FIG. 1. a) A schematic picture of the normal metal-superconductor junction in the 110 direction with a typical quasiclassical 
trajectory. 

b) A schematic graph of the pairing potential along the trajectory in a). 

Our approach was inspired by the study of a weak-coupling BCS superconductor modeled by a 
Fermi sea with an attractive local four-fermion interaction 

V^}(r)^(r)^(r)^ T (r); 

V < is the coupling constant. One decouples this interaction by the Hubbard-Stratonovich (HS) 
transformation, which introduces the pairing field A. The instability to superconductivity can then 
be detected by going to T = and looking at the energy change due to the opening of the gap |A|. 
The occupied single-particle states lower their energy as | A| 2 In | A| while the system raises its energy 
by |A| 2 /(— V), which is the extra term in the hamiltonian introduced by the HS transformation. 
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Since the non-analytic decrease wins over the analytic increase as |A| — > 0, an arbitrarily weak 
attractive V will, at low enough temperature, make the system superconducting. The decrease of 
the single-particle energies is non-analytic due to contribution from the states initially close to the 
Fermi surface. We conclude, therefore, that these low-energy states drive the BCS transition. 

A similar argument works for the d — > d + is transition. If there is an attractive interaction in the 
s-wave channel with strength V s < 0, we may use the HS transformation to introduce the s-wave 
component of the pairing field A s on top of the dominant A^, leading to an extra positive term in 
the hamiltonian | A s \ 2 /(— V s ), just as in the BCS case. When the "zs" component appears close to 
the surface, the shape of the order parameter along the trajectory in Fig. 1 changes, as shown in 
Fig. 2. To study the instability to the transition, we again look at T = and |A S | — ► 0. The energy 
of the ZBSs []to the lowest, that is first, order in A s then changes to fll5 | 



-oo 



E e [A s } = ± / d P 2\f(9,p)\ 2 lmA s , (1) 



-oo 



where the upper (lower) sign corresponds to the up- (down-)moving trajectory. Here 

.9(0, P). 

is a solution of the Andreev equation along the trajectory in the direction 9 at the point p. For the 
ZBSs, g(9, p) = =R/(0) p)- At T = 0, only the down-moving ZBSs will be occupied, and their energy 
will decrease linearly. It turns out |15[ that the energy changes of the remaining states (non-ZBSs) 
cancel each other out. Thus, upon the appearance of small |A S |, the energy of the occupied (zero- 
bias) states is lowered linearly, which, for small enough |A S |, wins over the quadratic increase of the 
HS term for arbitrarily weak s-wave attraction V s . It follows that the transition d — > d + is is driven 
by the ZBSs. 



1 We shall abuse the terminology and call these states ZBSs even after their energy has been shifted away from zero. 
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FIG. 2. The pairing potential along the trajectory in Fig. la). The corresponding twist of the phase of the order parameter 
is clockwise. 

For YBCO, the experiments give |A^| ~ 30 meV (amplitude) and |A S | ~ 1 meV, so we assume 
the first-order perturbation formula (|T|) holds up to the experimental value of |A S |. We can then 
calculate |A S | by minimizing the energy of the system. When we sum over all occupied ZBSs, we 
find that the energy per CuO plane per unit length of the surface is 

E[s(x)} = Jdx^r+ I ^-d9cos9E e [s(x)}, (2) 
where s(x) = ImA s (x). Since s extends into the bulk only as far as the <i-wave coherence length 



£ = Hvp/\Ad\, we can estimate 



E[s] ~ -T^rt ' k F s, (3) 



which gives 

(4) 

The experimental values s ~ 1 meV, kp ~ 1 A -1 , and £ ~ 10 A give us \V S \ ~ 10 meV A 2 . (Since 
the CuO plane is two-dimensional, V s has dimension EL 2 rather than EL 3 .) 



The full variational calculation is presented in We only remark here that the solution of the 



variational equation obtained from (^) agrees with the contribution to A from the occupied ZBSs in 
the gap equation 



° k 



A ZBS (x) = i(-V s ) J ^d92\f(6,x/cos6)\ 2 , (5) 



-tt/2 



that is, Azbs(%) = is(x), which shows the internal consistency of the picture. Unlike the BCS gap 
equation, @ is an explicit formula for Azbs(x); there is no Azbs on the right-hand side. The 
physical reason for this is that within the first-order perturbation theory, it is only the sign of s(x) 
(and not its magnitude) that determines which ZBSs are occupied. 

Our argument has shown only that d + is is the favorable state at T = 0. To see what happens at 
finite temperatures, we need to minimize the free energy F[s(x)], which is obtained from (§) when 
we replace E 9 [s(x)} by (— T) ln(l + exp(— E e [s(x)]/T)), making the lowest order quadratic in s(x). 
The order-of-magnitude estimate gives 

F M~(rb-T) s2+0 < s4 >- < 6 > 

Hence, the system is unstable to the transition to the d + is state even at finite temperatures. The 
transition temperature is 

T, ~ — J— , (7) 

which is of the same order of magnitude as |A S | T=0 (see (f|)). 

Now we shall go back to T = and observe that the presence of the is component induces a twist 
in the phase (p of the order parameter from —n to for an up-moving trajectory and from to 7r for 
a down-moving one, see Fig. 2. This twist implies a current flowing down; its density at the point 
p of a quasiclassical trajectory labeled by 9 is equal to 

j(9,p) = ^-n^dM0,p), (8) 

where = kp/rc is the density of the one-dimensional Fermi sea {kp is the Fermi wave vector). 
Current flowing down is also expected from our previous argument that showed that only the down- 
moving ZBSs are occupied. Each of them contributes 

j(9,p) = ev F (\mp)\ 2 + \9(e,p)\ 2 ) (9) 

to the current density at a given point; it turns out again |15| that the contributions from the non- 
ZBSs cancel. In both approaches, we get the total current by summing up contributions from all 
the quasiclassical trajectories. The result is 
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I = ^L (10) 

per CuO plane in both calculations. To get an order-of-magnitude estimate, we take vf ~ 10 5 m/s 
and kp ~ 10 10 m -1 , so I ~ 10~ 5 A. The agreement shows that all the current is carried by the ZBSs. 
It is important that not only the magnitude, but also the direction of the current, agrees in both 
calculations. The reason for this agreement is that the ZBSs moving in the direction of the current 
are shifted down in energy and thus occupied, whereas those moving against the current are shifted 
up and unoccupied. We wish to stress that this is exactly opposite to the sign of the Doppler shift: 
the states moving along the current would be Doppler-shifted up, whereas those moving against the 
current would be Doppler-shifted down. 

In conclusion, the first-order perturbation theory shows that the d — * d + is transition is driven 
by the ZBSs and that the transition temperature is finite and of the order |A s |t=o- The ZBSs that 
are being pushed down by the is component carry surface current; they are not Doppler-shifted by 
this current. 

I want to give special thanks to M. Stone for suggesting the problem and discussing it with me. 
I have also benefited from discussions with H. Aubin, L. H. Greene, A. J. Leggett, D. E. Pugel, R. 
Ramazashvili, M. Turlakov, and H. Westfahl. I would like to thank C. M. Elliott for proofreading 
the manuscript. The project was supported by the grant NSF-DMR-98-17941. 
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